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We introduce a gauge invariant definition of a monopole on the lattice. Tlie construction is based on the 
observation that for each Wilson loop there exists an extra U{1) group which leaves the loop invariant. Since the 
lattice formulation utilizes the language of Wilson loops, the definition of the monopole charge in terms of this 
plaquette dependent U{1) is gauge invariant. The explicit construction of gauge invariant monopoles is presented 
both in continuum and on the lattice. 



Introduction of magnetic monopoles in non- 
Abelian pure gauge models has a long history (see 
jl] and references therein). The basic problem 
is that the monopoles are intrinsically U{1) ob- 
jects and, in absence of the physical Higgs field, 
the choice of the U{1) turns out to be a mat- 
ter of gauge fixing. The standard procedure^ to 
define monopole charge is to partially (up to a 
remaining U{1)) fix the gauge and then construct 
Abelian monopoles in this particular gauge. Due 
to the ambiguity of gauge fixing prescription the 
construction is not unique and it is a separate 
question which monopole definition is physically 
relevant. 

1. The crucial observation |^ which is used 
in considerations below is that in lattice formu- 
lation there is in fact a natural C/(l) subgroup 
locally embedded into SU{2). Indeed, the action 
of SU{2) lattice gluodynamics is constructed in 
terms of elementary Wilson loops Up = e*^p°'/^ — 
^i\Fp\npcr/2 ^ Tip = I and therefore possesses an ad- 
ditional symmetry Fp —^ Fp + Airfip which is evi- 
dently lost in the naive continuum limit. In other 
words Wilson loop itself defines a natural U{1) as- 
sociated with it as the group of rotations around 
ftp. Since these Abelian rotations are defined in 
terms of a Wilson loop, the definition of the U (1) 
is gauge invariant by construction. The appli- 
cation of this procedure to each plaquette gives, 
basically, a gauge invariant fixation of t/(l) which 
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is plaquette dependent. 

To make the next step and define monopoles 
in terms of this local, gauge invariant U{1) we 
utilize the interpretation of Wilson loop as an 
evolution operator of quantum mechanical system 
1^. As follows from the standard representation 
W{C) = W{T) = P exp{i /J Adt} the Hamilto- 
nian governing evolution is if = —A and quan- 
tum mechanical state space coincides with irre- 
ducible representation space of SU{2). Here C 
is an arbitrary closed contour parameterized by 
t £ [0;T]. Furthermore, we introduce the con- 
ventional [Q coherent state basis parame- 
terized by a set of unit three-dimensional vectors 
n, "n? = \. The action of an arbitrary element 
g e SU{2) in this basis ii^ g \n) ^ e"^ \ n' ). A 
special case of cyclic evolution 

W{T) I n(0) ) = e*'^(^) | n{0) ) (1) 

is of particular importance. In Eq. (|^) the evolv- 
ing state I n{t) ) is determined by 

P exp{i / Adr} I n(0) } = e"^^'^ \ n(t) ) . (2) 
Jo 

Clearly the cyclic state | n ) always exists and 
moreover ^TiW(T) = cos(p{T). From (||) it di- 
rectly follows that 

^iT) = £ (^{n\A\n.)+^{n\^^\n)yt= (3) 

^ For simplicity we consider fundamental representation 
only, generalization to higher representations is straight- 
forward. 
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where the state | n ) has been smoothly extended 
from the contour C into an arbitrary surface Sc 
bounded by C. Note that Eq. (^) cannot be used 
in the Wilson loop calculation since construction 
of the evolving state | n ) requires knowledge of 
the Wilson loop itself. Nevertheless, Eq. (H) is 
useful theoretically since it represents the phase 
angle of Wilson loop as an integral of abelian 
't Hooft tensor Note also that the angle ip{T) 
is well defined only modulo 27r. Indeed, one can 
verify explicitly that is changing — > + 27rfc, 
k £ Z under gauge transformations. From now 
on we always take ip(T) in "fundamental domain" 

— TT < if < TT. 

Consider now the infinitesimal version of 
Eq. (||) when contour C is the boundary of el- 
ementary surface element Sa 



(4) 



(p ~ (j) A (fiA) +n Adn A dnj Sa . 

It is straightforward now to integrate 
bitrary closed two-dimensional surface S'^/^y^ in 
physical space. Usual assumption of fields conti- 
nuity guarantees that r? is a smooth field on S^i^^^ . 
Therefore integral is non-zero in general 



over ar- 
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if = 2t:Q . 



(5) 



due to the second term in Eq. (|j) which is widely 
known as Berry phase In fact Eq. (||) calcu- 
lates the gauge invariant monopole charge Q con- 
tained inside S^f^y^. Gauge invariance is evident 
since integral is constructed in terms of infinites- 
imal Wilson loops only. The term "monopole 
charge" is also quite natural because Eq. (||) co- 
incides with well known Abelian monopole con- 
struction when gauge fields have the same color 
structure. Moreover, Eq. being considered 
in four-dimensional space-time defines a closed 
world-lines of topological defects which have a 
natural interpretation of monopole trajectories. 

To conclude this section we would like to em- 
phasize that gauge invariant monopole charge 
definition is only possible when gluodynamics is 
considered as limiting case of lattice gauge mod- 
els because of periodicity of gauge fields action 



on the lattice (or more details see One 
can verify that Eq. (||) considered in the context 
of conventional continuum gluodynamics is trivial 
and produces identically zero magnetic charge. 

2. Being extremely simple in continuum limit, 
Eq. is quite non-trivial to implement on the 
coarse lattice. The problem is that lattice dis- 
cretization is not suitable to calculate most of 
topological invariants (e.g. instanton number). 
In our case the problem is even more severe 
since effective Higgs field in Eq. (||j5|) is de- 
fined not in lattice sites as usual, but on ele- 
mentary 2-cells (plaquettes). Below we illustrate 
the monopole charge calculation in case of single 
three-dimensional cube. More realistic calcula- 
tions will be presented elsewhere [Q . 

To begin with we consider a single plaquette 
situated at lattice site x and directed along /z, 
V space-time directions. It is straightforward to 
calculate plaquette matrix Up as ordered product 
of links. For fundamental representation there 
are two eigenvectors | n± ): Up\n±) = e^^'^\ n± ) 
which are related by n+ = — n_. While only 
a single plaquette is considered there is no way 
to distinguish between | n_|_ ) , | ri_ ) and one can 
take either of them as initial state \ ni) to be as- 
cribed to point 1 on the plaquette (see Fig. 1). 
Starting from | ni ) one builds the corresponding 
evolving states | ) in all other plaquette cor- 
ners: Up_{x) \ fii ) = e*'''^ I n2 ) etc. Thus the lattice 
implementation of Eq. (||,||) is obtained. 
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Fig. 1 Fig. 2 

Consider now the intersection of two plaque- 
ttes. Fig. 2. Evidently, on coarse lattice the states 
build separately on each plaquette differ drasti- 
cally. In particular, the states \ni), i = 1,2 
and i = 3, 4 nearest to common link do not form 
continuous vector field ft which is needed to go 
from Eq. (^ to Eq. (||). Moreover, there is also 
a mentioned ambiguity in the choice of initial 
states I ft± ) . Since the problem is only due to 
the lattice coarseness there should be no differ- 
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ence in continuum limit between various ways to 
overcome it. We propose to introduce additional 
two dimensional cells in the intersection of ev- 
ery two plaquettes which do not lie in the same 
plane, Fig. 3. The newly introduced hnks V12, 
V34 may be defined in fact unambiguously. In- 
deed, among various SU{2) matrices which sat- 
isfy V12 1 n2 ) = e*" I ni ) there is only one which 
corresponds to geodesic motion n2 —>■ ni. Note 
that these new links are in adjoint representation 
since for example fii and 722 transform in the same 
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way under gauge transformation: | ) — > g\ni), 
i = 1,2. The ambiguity in choice of initial state 
on every plaquette may be resolved in similar 
manner. Namely, one can impose a requirement 
that resulting field n should be the smoothest one. 

To test the proposed algorithm we have consid- 
ered the simplest case of single three-dimensional 
cube. The first test was to generate random 
pure Abelian gauge fields, calculate the monopole 
charge in standard way and then apply our 
method after random non-Abelian gauge trans- 
formation. We found that our method works per- 
fectly on single 3-cube reproducing the known 
monopole charge in all cases. Furthermore, we 
have found a reasonable behavior when non- 
Abelian gauge fields Ui = const - {1 + CiSf^QeiCT*} 
were generated (here Si G are random 

numbers). On Fig. 4 the average monopole charge 
as a function of C is shown. 

3. We have proposed a gauge invariant defini- 
tion of Abelian monopole charge in gluodynamics 
which is turn is only possible when gluodynam- 
ics is considered as limiting case of lattice gauge 
models. The crucial element of the presented con- 
struction is the observation that actual symmetry 
group of SU{2) LGT is in fact SU{2) x U{1) and 




therefore monopoles defined with respect to sec- 
ond factor are SU{2) invariant. The explicit cal- 
culations are carried out both in continuum and 
on the lattice. Being transparent in continuum 
limit the actual calculations become rather intri- 
cate on the coarse lattice. We have proposed a 
way to overcome this difficulty which allows to 
investigate the monopole dynamics numerically. 
Unfortunately, we were not yet able to implement 
the method in realistic calculations since it re- 
quires changing of lattice geometry and standard 
Monte-Calro algorithms. The work in these di- 
rections is currently in progress. 
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